Bogoliubov- Cerenkov radiation in a Bose-Einstein condensate 
flowing against an obstacle 
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We study the density modulation that appears in a Bose-Einstein condensate flowing with superso- 
nic velocity against an obstacle. The experimental density profiles observed at JILA are reproduced 
by a numerical integration of the Gross-Pitaevskii equation and then interpreted in terms of Cerenkov 
emission of Bogoliubov excitations by the defect. The phonon and the single-particle regions of the 
Bogoliubov spectrum are respectively responsible for a conical wavefront and a fan-shaped series of 
precursors. 
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The Cerenkov effect was first discovered in the electro- 
magnetic radiation emitted by charged particles travel- 
ing through a dielectric medium at a speed larger than 
the medium's phase velocity [1]. A charge moving at 
the speed v is in fact able to resonantly excite those 
modes of the electromagnetic field which satisfy the kine- 
matic Cerenkov resonance condition uj era {k) = v • k: 
part of the kinetic energy of the particle is then emit- 
ted as Cerenkov radiation, with a peculiar frequency and 
angular spectrum [2[. Electromagnetic waves in a non- 
dispersive medium of refractive index n have a linear dis- 
persion law relation uo em (k) — ck/n: the Cerenkov con- 
dition is then satisfied on a conical surface in /c-space of 
aperture cos</> = c/(nv), which corresponds to a conical 
wavefront of aperture = tt/2 — <p behind the particle. 
Thanks to the interplay of interference and propagation, 
much richer features appear in the spatial and /c-space 
pattern of Cerenkov radiation in dispersive media 0, Q 
and photonic crystals [5j]. 

The concept of Cerenkov radiation can be general- 
ized to any system where a source is uniformly moving 
through a homogeneous medium at a speed larger than 
the phase velocity of some elementary excitation to which 
the source couples. Many systems have been investigated 
in this perspective, ranging from e.m. waves emitted by 
the localized nonlinear polarization induced by a strong 
light pulse travelling in a nonlinear medium 0, Q , to the 
sonic waves generated by an airplane moving at super- 
sonic velocities, to phonons in a polaritonic superfluid [8[, 
and in a broader sense, to the surface waves emitted by a 
boat moving on the quiet surface of a lake Q . In this Let- 
ter we compare our theoretical results of the density per- 
turbation induced in a Bose-Einstein condensate (BEC) 
which flows against a localized obstacle at rest with the 
experimental images taken by the JILA group • Mod- 
ulo a Galilean transformation, the physics of a moving 
source in a stationary medium is in fact equivalent to the 
one of a uniformly moving medium interacting with a sta- 
tionary defect. The experiment has been performed by 
letting a BEC expand at hypersonic speed against the lo- 



calized optical potential of a far-detuned laser beam. The 
observed density profiles are successfully reproduced by 
numerically solving the time-dependent Gross-Pitaevskii 
equation and physically interpreted by a simple model of 
Cerenkov emission of Bogoliubov excitations by a weak 
defect. 





FIG. 1: Experimental [T(| density profiles (integrated along z) 
of a BEC hitting an obstacle at supersonic velocities v/c s = 13 
(a) and 24 (b). The angles of the conical wavefronts are 
sin(0) = 0.73 and sin(0) — 0.43, respectively. The conden- 
sate flow is from the right to the left. 



The experiment. The experimental results analyzed in 
this paper have been obtained by the JILA group [To[ 
with a gas of N = 3 x 10 6 Bose-Einstein condensed 87 Rb 
atoms confined in a cylindrically symmetric harmonic 
trap of frequencies {cu r ,cu z } = 27r{8.3, 5.3} Hz. The BEC 
is slightly cigar shaped, with the long axis pointing in 
the direction of gravity (z-axis), and a Thomas- Fermi 
radius in the x — y plane of 31.6 /mi. At t = 0, the sign 
of the confining potential in the xy plane is inverted, so 
that the BEC undergoes an antitrapped expansion in the 
xy plane. At the same time, the confinement along z is 
switched off and replaced by a linear magnetic field gra- 
dient to cancel the effect of gravity. The obstacle consists 
of the repulsive potential of a blue-detuned dipole beam 
with a round Gaussian profile. The beam is directed in 
the z-direction and placed in the vicinity of the trapped 
condensate, and remains in this position during the ex- 
periment. Images of the BEC density profile after differ- 
ent expansion times t exp are then taken by means of de- 
structive absorption imaging. Two examples are shown 
in Fig.(pQ). The field of view is centered in the region 
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around the defect in order to observe the perturbation 
that it creates in the expanding condensate. In this re- 
gion, the condensate is flowing from the right to the left. 
During the expansion, the condensate is strongly acceler- 
ated by the inverted harmonic trap, so that the local den- 
sity at the defect position decreases in time, and the local 
speed v increases. As a result, images taken after increas- 
ing expansion times correspond to increasing values of the 
Mach number v/c s , c s being the local value of the sound 
velocity at the defect position. In FigQJa,b), the expan- 
sion times are t = 33.6, 50 ms, long enough for the flow 
to be strongly supersonic v/c s = 13, 24. Downstream 
of the defect the atomic density is strongly reduced in a 
conical shadow region, separated from the unperturbed 
condensate region by a conical wavefront followed by a 
fan-shaped series of precursors which extends far in the 
upstream direction. For increasing Mach number v/c 8 , 
the aperture of the shadow area decreases, as well as the 
spatial period of the precursors. Physical interpretations 
for all these features will be provided in the following. 

Gross-Pitaevskii simulations. As a first step, the ex- 
perimental data have been reproduced by a numerical 
simulation of the 3-dimensional, time-dependent Gross- 
Pitaevskii equation by means of a real-space-product fi- 
nite element discrete- variable-representation [l]J, a nu- 
merical method that has already seen effective service 
in solving several large-scale computational problems in- 
volving BECs [l2[. Results for the exact parameters of 
the experiment are shown in Figj2^a,b). 

The laser obstacle consists of a repulsive potential with 
a Gaussian shape V(r) = Vmax exp(— 2r 2 /w 2 ). The max- 
imum height Vm ax = 4961 hw r of the potential is much 
larger than the chemical potential of the BEC in the ini- 
tial trapped configuration and its width w — 17 fim is 
~ 50 times larger than the BEC healing length: the ob- 
stacle can therefore considered, for all practical purposes, 
as a rigid cylinder. In Figs|2](a,b), the defect is situated 
at a distance of respectively d = 75, 87 /im from the con- 
densate center and the imaging times are t = 33.6, 50 ms. 
The density profiles emerging from the simulations are 
in agreement with the experimental findings of Fig.(pQ). 
The Mach numbers in the GPE simulations are 12.4 (a) 
and 23.8 (b), and the angles of the conical wavefronts 
delimiting the shadow areas are sin(0) = 0.7 (a) and 
sin(0) = 0.48 (b). The values of the angles are in excel- 
lent agreement with the experimental ones: sin(0) = 0.73 
and sin(0) = 0.43 for a Mach number 13 and 24, see 
Fig.(pQ). The GPE simulations also reproduce the fan 
of precursors moving upstream of the defect [13]. By 
comparing panels (a) and (b) of FigEl one sees that the 
longer the imaging time, the longer the spatial distance 
the precursors have traveled. 

Both the spatial period of the precursors, and the aper- 
ture of the cone decrease for an increasing Mach number 
v/c 3 , but this latter is always significantly larger than 
the Mach angle such that sin(0) = v/c s . As the flow ve- 



locity of the expanding BEC is not uniform, but rather 
radial, the geometric angle formed by the tangent of the 
obstacle through the BEC center and the line joining the 
obstacle and the BEC centers has to be added to the 
Mach angle. In addition to this, the geometrical shadow 
effect provides a simple, ballistic, explanation also for the 
strong density reduction inside the cone. This interpre- 
tation has been verified by repeating the simulations for 
obstacles of smaller size: as it is shown in Figged), the 
aperture of the cone significantly decreases and asymp- 
totically approaches the Mach angle in the limit of small 
defect, see Fig. ((3]). 
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FIG. 2: GPE density profiles integrated along the ^-direction 
obtained with different expansion times and sizes of the ob- 
stacle : a) t = 33.6ms, w=17 /mi.; b) t = 50ms, w=17 /im; c) 
t 38.4ms, w=3.74 /mi; d) t 38.4ms, w=0.374 /mi. 



FIG. 3: Conical wave front aperture as a function of the ob- 
stacle width (in "trap units", 1 t.u. =3.74 /am) calculated 
with GPE. The green line is the Mach angle. 



Bogoliubov-Cerenkov theory. A simple way of getting 
a qualitative physical understanding of the experimen- 
tal and numerical results is to use an approximate the- 
ory where the condensate is assumed to be homogeneous 
and uniformly flowing at a constant speed vo and the 
defect, localized around x = 0, induces a weak density 
perturbation 0, [TH, [l6[ • Under this approximations, an- 
alytical results for the stationary solution of the GPE can 
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be obtained by means of the Bogoliubov linear response 
theory. As boundary condition, an incident plane wave 
of the form ^o( x ? t) = ipo e *( k o*-^o£) i s taken, where we 
have set foko = mvo, po = \^o\ 2 , = + gpQ. 

Here g = A-K?i 2 a s /m in terms of the atomic s-wave scat- 
tering length a s and mass m. The sound velocity is 
c s y/gpo/m. 

The BEC wavefunction ^(x, t) = [^o(x) + 
is the sum of the unperturbed 
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*V(x,t)] 

wavefunction iJjq (x, t) plus the small perturba- 
tion 5ip(x.,t). Introducing the compact notation 
6ip(x.,t) = (J^(x, t), ^*(x, t)) T , the evolution can be 
written in the form of a linear equation 
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with a source Fd proportional to the defect potential V d : 



F d (y) = V d (y) 
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The Bogoliubov matrix C has the usual form: 
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and its eigenvalues give the energies of the elementary 
excitations of the system around the unperturbed sta- 
tionary state. In our specific case of a homogeneous and 
homogeneously flowing BEC, the Bogoliubov dispersion 
law is 



cj(k) = v-(k-k )± 



ft(k-ko) 2 /ft(k-ko) 2 
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(4) 

As we can see in the dispersion spectra shown in the left 
column of FigHl the main effect of the flow consists of 
the additional term vo • (k — ko) which tilts the dispersion 
and add vo to the propagation group velocity of all the 
Bogoliubov excitations. The ± branches correspond to 
respectively the particle- and the hole-like branches of 
the Bogoliubov dispersion, and are images of each other 
under the transformation k — » 2ko — k, co — > -co. The 
steady state in the presence of the defect potential V d is 
obtained from the motion equation (pQ) as: 



5$ d = -(C-iO + l) 
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The infinitesimal imaginary term is required to satisfy the 
boundary condition that at t = — oo no Bogoliubov exci- 
tations were present. As can be easily seen from the reso- 
nant denominator, the time- independent defect potential 
excites the Bogoliubov modes whose energy is o;(k) = 0. 
Graphically, these k modes can be identified in the plots 
in the left column of FigH] as the intersection points of 
the dispersion surface u(k) with the uo = plane. Mod- 
ulo a Galilean transformation, this condition corresponds 



to the usual Cerenkov resonance condition [2]. Depend- 
ing on whether the flow speed is slower or faster than the 
speed of sound in the BEC, two regimes can be identified. 

In the subsonic regime < c s , no intersection exists at 
k 7^ ko, which physically means that no Bogoliubov mode 
can be resonantly excited by the defect. The BEC is su- 
perfluid and can flow around the defect without suffering 
any dissipation, in agreement with the Landau criterion 
of superfluidity. 
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FIG. 4: Dispersion of Bogoliubov modes (left column), curve 
T in k space describing the resonantly excited Bogoliubov 
modes (center column) and spatial density profile (right 
colomn). The supersonic flows are vo/c s = 1.1, 1.5, 3 in each 
row, respectively. 



On the other hand, in the supersonic regime, vq > c Sl 
the set of k vectors satisfying o;(k) = is not empty, 
but rather corresponds to the closed curve Y shown in 
the central column of FigHJ Some of the kinetic energy 
associated to the flow is therefore dissipated as radiation 
of Bogoliubov modes. In real space, the perturbation ra- 
dially propagates from the defect with a velocity fixed 
by the group velocity of the mode v g = Vk^. For each 
mode k, the direction of v g corresponds to the outward 
normal to the curve T. Since we are here considering a 
stationary state, the density perturbation propagates to 
infinity in all available directions. Yet, even though T is a 
closed curve, the possible directions of the group velocity 
do not fill the whole unit sphere: because of the singu- 
larity at ko there exists a region of excluded directions, 
which in the plots in the left column of FigH] correspond 
to the unperturbed region inside the Cerenkov cone. 

The spatial profiles of the density perturbation pre- 
dicted by the linear theory can be understood by an- 
alyzing the different regions of T. Close to the singu- 
lar point at k = ko, the curve consists of two straight 
lines separated by an angle 2<p defined by the usual 
Cerenkov condition cos(/> = c s /vq. This region corre- 
sponds to the small wavevector part of the Bogoliubov 
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spectrum, which is indeed linear and allows for an in- 
terpretation of the corresponding radiation in terms of 
the usual Cerenkov effect in non-dispersive media. In 
the plots in the right column of figlH the strongest den- 
sity perturbation lies on the Mach cone of aperture 9 
such that sin# = c s /vq, and directed in the downstream 
direction; being the flow velocity larger than the sound 
velocity, the density perturbation is dragged away from 
the defect. The rounded region of T for k nearly anti- 
parallel to ko is absent in the standard theory of the 
Cerenkov effect in non-dispersive media, and originates 
from the quadratic, single-particle-like, part of the Bo 
goliubov spectrum. The group velocity at the right-most 
point of T is equal to v g = \c 2 s /vq — vq\ so that the pertur- 
bation propagates in the upstream direction with respect 
to the BEC flow and gives an oscillating density modula- 
tion of wave- vector |k — ko| = 2m(vQ — c 2 ) 1 ^ 2 /h. This is a 
direct consequence of matter-wave interference between 
the incident BEC and the scattered wave off the obsta- 
cle [ItJ- In the intermediate points on the curve, there 
is a one-to-one correspondence between excited Bogoli- 
ubov modes of wavevector k and the direction x such 
that x is parallel to the group velocity v p (k). For each 
point x = xx in space, this correspondence determines 
the wavevector k of the Bogoliubov mode that is able to 
reach it, and consequently the wavevector k m = ko — k of 
the density modulation in its vicinity. Despite the under- 
lying approximations, the features of the Bogoliubov the- 
ory qualitatively reproduce the full GPE simulations and 
therefore clarify the interpretation in terms of Cerenkov 
waves of the experimental findings. 

In conclusion, we have presented a combined theoret- 
ical and experimental study of the response of a Bose- 
Einstein condensate flowing against a localized obsta- 
cle at supersonic velocity [18|- A numerical study of 
the full Gross-Pitaevskii equation has been performed 
to fully simulate the experiment. The findings are in 
agreement with the experiment, and have been inter- 
preted within the Bogoliubov approximation in terms of 
the Cerenkov emission of phonons by the defect: this cre- 
ates a Mach cone density modulation downstream of the 
defect, as well as an additional upstream fan-shape per- 
turbation, which originates from the interference between 
the particle-like Bogoliubov excitations and the underly- 
ing BEC. For growing strengths and sizes of the obstacle, 
the numerical study shows a continuous evolution of the 
density modulation towards a soliton train [19[ , as typical 
of the dynamical evolution of shock- waves in dispersive, 
dissipationless nonlinear ID wave equations [2o| . 
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